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Abstract. In the present paper we extend the method to detect Pomeranchuk 
instabilities in lattice systems developed in previous works to study more general 
situations. The main result presented here is the extension of the method to include 
finite temperature effects, which allows to compute critical temperatures as a function 
of interaction strengths and density of carriers. Furthermore, it can be applied to 
multiband problems which would be relevant to study systems with spin/color degrees 
of freedom. Altogether, the present extended version provides a potentially powerful 
technique to investigate microscopic realistic models relevant to e.g. the Fermi liquid 
to nematic transition extensively studied in connection with different materials such 
as cuprates, ruthenates, etc. 



PACS numbers: 71.10.Hf, 71.10.Fd, 71.10.Ay 



Pomeranchuk instabilities in lattice systems 



2 



1. Introduction 

The question of the stability of the two-dimensional Landau Fermi liquid is long-standing 
and of fundamental importance. One important motivation is the non-Fermi liquid 
behavior observed in high-Tc superconductors above the critical temperature, which 
has in part promoted in the last years a surge of interest in possible Pomeranchuk 
instabilities pQ|2i[3j|IJ[5j[6j[7J[8j|9j [10] . In particular, the nematic transition due to a 
Pomeranchuk instability in two dimensions has been studied in detail in [U [UJ [121 02] ■ 
For a recent review with an extensive list of references see [H]. Later on, similar 
ideas were applied in [151 EH EL7] to include spin degrees of freedom. In connection 
to ruthenates [T8l \W\ |20"1 |2T] . Pomeranchuk instabilities have been advocated in order 
to explain experimental observations in the metamagnetic bilayer ruthenate Sr 3 Ru207, 
where an intermediate nematic phase is observed in a magnetic field. Another system 
where Pomeranchuk instabilities may play a role is in the hidden order transition 
observed in URu2Si2 [22], in this case in connection with spin-antisymmetric interactions. 

Motivated by these observations and ideas, we have performed a systematic study 
on how to analyze Pomeranchuk instabilities in lattice systems, with the aim to account 
for the lattice effects, which are expected to be strong in most of the strongly correlated 
electron systems 0. The effects of the lattice on the analysis of Fermi liquid instabilities 
have been previously considered in e.g. P HI [71 El [131 [231 Ell EHl [261 Ell 

Instead of parameterizing Fermi surface (FS) deformations using a basis of functions 
of the point group symmetry of the lattice, which would be case dependent and 
cumbersome in general, in previous works we developed a method that is formally 
independent of the underlying lattice and straightforward to be applied [301 EI]- This 
method can be henceforth used to detect Pomeranchuk instabilities in a generic model 
defined in an arbitrary lattice, which for completeness is reviewed in Section [2j 

In order to get closer to experiments, we extend the method to include finite 
temperature effects. This allows to study the extended phase diagram and in particular, 
to compute critical temperatures as a function of the chemical potential and interaction 
strengths. On the other hand, in order to study models relevant to the systems 
mentioned above, we need to include extra structure such as spin or band index, which we 
show how to handle by identifying the additional degrees of freedom that parameterize 
spin/color flipping excitations. 

The paper is organized as follows: in Section [2] we review the method developed 
in [301 EI] • Then in Section E] we present the details of the extension of the method 
to include finite temperature effects. In Section [4] we show an example of application 
of the extended method to the finite temperature analysis of a model with an s-wave 
instability. In Section [5] we show that internal degrees of freedom, such as spin or band 
indexes, become manifest as an extra matrix structure on the space of deformations, 
which has to be incorporated in the analysis of the instabilities. 

| Except for the case of the two-dimensional electron gas in strong magnetic fields 
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2. Brief review of the spinless case at zero temperature 

According to Landau's theory of the Fermi liquid, the change in the Landau free-energy 
Q = E — TS — fiN as a functional of the change 5n(k) in the equilibrium distribution 
function at finite chemical potential \x can be written, to first order in the interaction, 

as 



Here e(k) is the dispersion relation that controls the free dynamics of the system, the 
interaction function /(k, k') can be related to the low energy limit of the two particle 
vertex. Note that we are omitting spin indices, and considering only variations of the 
total number of particles n(k) = n-j-(k) + ^(k). This implies that the considerations 
that follow in the present section will be valid in the absence of an external magnetic 
field and at constant total magnetization. 

Starting from this description, Pomeranchuk has shown long time ago a simple way 
to detect instabilities of the Landau Fermi liquid theory by considering the change in sign 
of the energy when slightly modifying the occupation numbers [32] . In his original paper 
he worked with a three dimensional model with a spherical FS, and hence computations 
could be easily performed by using spherical harmonics. 

In [30| l3Tj we have extended the original Pomeranchuk idea to deal with lattice 
systems or, more generally, to systems with a Fermi surface with an arbitrary shape 
where no base functions, playing the role of the spherical harmonics in Pomeranchuk's 
case, can be easily identified. This problem has been previously studied using 
renormalization group techniques in [23] for the square lattice case with next-nearest- 
neighbor hopping, where it was reduced to the diagonalization of a matrix whose rank 
is of the order of the number of sites. 

In the present paper we generalize our approach to deal with finite temperature 
effects. We also include spin degrees of freedom as well as any other internal index, such 
as the band index in a multiband system. 

The key idea in Pomeranchuk's analysis is to characterize deformations 5n(k) that 
would lead to 5Q < 0, then corresponding to an instability. In order to perform this 
analysis in an arbitrary lattice, the first step is to change variables in (JTJ from (k x , k y ) 
to normal Fermi variables (g, s) defined as 



where the new variable g is normal to the FS, which is defined by g = 0, and the 
tangent variable s varies between — ir and ir. The associated Jacobian is given by 
J~ 1 (g, s) = \d(g, s)/dk\ . We will assume that this change of variables is well defined 
almost everywhere in momentum space, with the possible exception of isolated points 
at which the density of states may diverge (van Hove singularities). 

In a stable system, the energy (JTJ should be positive for all d~n(k) that satisfy the 




(1) 



s = s(k) , 



(2) 
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constraint imposed by the Luttinger theorem 

Jd 2 k5n(k) = 0. (3) 

Pomeranchuk's method roughly consists on exploring the space of solutions of the 
Luttinger constraint (jHD to find the possible 5n(k) that turn the energy (CQ) into negative 
values, thus pointing to an instability of the system. 

In the new variables, an arbitrary deformation can be parameterized by 5g(g, s) 
and we can write Sn(g, s) at zero temperature as 

8n(g,s) = H[g + 5g(g,s)}-H[g], (4) 

where H is the unit step function. Replacing this into the constraint (131), changing the 
variables of integration according to ([2]) and performing the integral to lowest order in 
Sg we get 



dsJ{s)5g{s) = 0. (5) 

Here J(s) = J(g,s)\ g=0 and Sg(s) = 8g(g, s)\ g=0 . In the case in which the FS has a 
nontrivial topology, the integral would include a sum over all different connected pieces. 

We see that in order to solve the constraint, we can write an arbitrary deformation 
Sg(s) as 

J(s)5g(s) = d s \(s) , (6) 

in terms of a free slowly varying function A(s). 

Using the change of variables (j2D and with the help of eqs. (jlj and (jHJ), we can 
rewrite the variation in the energy Sfl to lowest order in Sg as 

5n=J ds'J ds V>(s')~ (j-\s)6(s - s') + f(s, s')) , (7) 

where we call d s \(s) = if)(s) and f(s, s') = f(g, s; g r , s')\ g=g > =0 . 

Note that the variation in the energy Sfl has two terms, the first one contains the 
information about the form of the unperturbed FS via J~ 1 (s), while the second one 
encodes the specific form of the interaction in f(s,s'). There is a clear competition 
between the interaction function and the first term that only depends on the geometry 
of the unperturbed FS. 

We see in (J7|) that the stability condition is equivalent to the positive definiteness 
of a quadratic form 

5n = {i/>,if>). (8) 

for any function if> in the space of square-integrable functions defined on the Fermi 
surface Z^fFS]. This was the main achievement in our previous paper: we have managed 
to express the stability condition in a completely general form which is independent of 
the underlying geometry. 

The rest of the procedure is completely technical and it is described in detail in 
\30\ [31] . The main idea is to diagonalize this quadratic form by applying the Gram- 
Schmidt orthogonalization procedure to an arbitrary starting basis of Z^fFS], with the 
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pseudo scalar product defined by (j7j). This results in a new orthogonal basis, that may 
eventually contain functions with negative pseudonorms. In that case, the corresponding 
variation on the energy (J5J) will be negative and we say that we have found an unstable 
channel. 

This summarizes the generalized Pomeranchuk method that can be applied to 
any two dimensional model with arbitrary dispersion relation and interactions. Some 
additional comments are in order: First, in some cases the solution of the constraint 
that defines the variable s can be extremely involved, and then an alternative route 
should be followed as described in [3T]. Second, the minimal length that can be resolved 
in phase space is determined by the total number of sites of the lattice. In consequence, 
we do not need to consider deformations whose characteristic length is smaller that such 
minimal length. Then if the starting basis is ordered in decreasing order according to 
some characteristic length of the basis functions, we can limit our procedure to a finite 
number of them. 

It has to be stressed that in its present form, the method is not suitable to detect 
first order phase transitions of the kind already predicted for some two dimensional 
systems [29]. This is due to the fact that the energy is expanded up to second order in 
the variations of the occupation numbers. Nevertheless, since our method determines 
sufficient conditions for a phase transition to occur, additional first order instabilities 
can only enlarge the unstable region. 



3. Finite temperature case 

To generalize our procedure to finite temperature, the variations in the mean occupation 
numbers 5n(g, s) in (j3J) have to be replaced by the corresponding expressions at finite 
temperature, namely the Heaviside function H in (jlj) has to be replaced by the Fermi 
distribution F at temperature T 

F ^ = e -J B T + 1 ( 9 ) 

then 

5n(jg,s) = F[g+6g(g,3)]-F\g], (10) 

where as before Sg(g,s) is a small perturbation parameterizing the deformation of the 
FS. Replacing into the constraint (J3J), changing the variables of integration according to 
(J2J) we have 

J dgds J[g, s] (F[g + 5g(g, s)] - F[g]) = (11) 

and expanding to first order in 6g(g,s) we get 

Jdgjds Fjg]J(g, s)8g(g, s) = (12) 

where F t9 denotes derivative of F with the respect to g. This is the finite temperature 
version of the constraint fl5]), that has to be solved in order to obtain the independent 
degrees of freedom. 
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To obtain a solution, we note that deformations Sg(g, s) satisfying 
/ ds J(g,s)Sg(g,s) = 



(13) 



constitute a particular set of solutions of ( 1T2|) . In what follows, we will restrict our 
search for unstable deformations to that particular set, keeping in mind that more 
general solutions can only enlarge the unstable region of phase space. The solution of 
( 113"]) can be written as 



Now, replacing ( flOl) into the energy of a perturbation (|T]), expanding to second 
order in 5g (g, s), and using the solution ( fl4l) we get 



where ip(g, s) = d s X(g, s). Note that the term linear in d s \(g, s) vanishes by making use 
of f fTB"j) . This expression is the finite temperature version of ([7]), and allows us to write 
the energy as a temperature dependent bilinear form 



The main difference with the zero temperature case is that now both normal Fermi 
variables g and s appear in the integrals and henceforth the rest of the procedure, 
i.e. the search for negative eigenvalues to diagnose an instability, has to be modified 
accordingly. More precisely, one has to write the energy as the bilinear form ( fl5l) using 
the full expressions for the Jacobian J(g, s) and the interaction f(g, s; g', s'). Since now 
the deformations ip(g,s) are functions defined on the whole momentum plane, one has 
to choose an arbitrary basis {ip n } n( zN now of the space of functions L 2 [R 2 }. 

The rest of the procedure remains unchanged, and the instabilities can be diagnosed 
just as before, by applying the Gram-Schmidt orthogonalization procedure to the chosen 
starting basis of functions to find the modes which have a negative pseudonorm. 

Even if at this point the method can be said to be complete, in practice the 
calculations needed may be very involved. As already mentioned in section [3j it may be 
very difficult to obtain the explicit form of the Jacobian for a given problem. Moreover, 
even if the Jacobian is known, the integrals involved in eq. ( fl5l) may be very time- 
consuming for a useful numerical implementation. On the other hand, it may not be 
clear whether an instability found at finite temperature survives in the zero temperature 
limit. In order to shed light on these issues, we will go further into the details of the 
present analysis. 

First, note that even if the functions F >g [g] and F ;gg [g] have noncompact support, 
they are strongly suppressed outside a region of width k^T around the origin. This 
implies that at low enough temperatures, the main contribution to the integral ( fl5l) 
comes from the values taken by the functions J(g, s), f(g, s; g', s') and ip(g, s) close to 
g = 0. In particular, since the maximum of the function F >g \g] sits at g — 0, it is 



$g(g,s) = J (g,s)d s X(s,g) 



(14) 




(15) 



(16) 
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natural to expand the remaining factors in powers of g. The interaction function can 
be expanded as 

f(g,s;g\s') = J2fnm(s,s')g n g ,m (17) 

n,m 

while the Jacobian takes the form 

J-\g,s)=Y,Jn 1 (s)g n (18) 

n 

Decomposing the deformation ip(g, s) as 

^(g,s) = J2Ms)g n (19) 



5tt = £ I dsds'Ms)\Y. ((n + m + r + l)J^(s)F n g +m+r S rt 5(s - s') 

n,m^ £ r ,t 



we get 



+ f rt (s,s')F^ r F^ +t j^') (20) 
where we called F n g the n-th momentum of the function F }3 [g], namely 

K = J d 9 FM9 n (21) 

It has to be noted that, since F t3 [g\ is an even function, only even momenta contribute 
to the sum above. Using this fact, the above expression can be rewritten as 

6n=Y,[ dsds'Ms)l E((2^ + l)J^ n - m {s)Ff5 jk 5{s - s') 

n,m J * j t k 

+ f2k-n,2j-m(s, s')F 2 g k F 2 A ^ m { S ') (22) 



Since F 2 g k scales like (fcsT) 2fc , at low enough temperature we can keep only the lowest 
orders in k and j in the formula. On the other hand, since J~ l (g,s) and f(g,s;g',s') 
are assumed to be regular at g = 0, they contain only non-negative powers of g. Then 
in the above sum 2k — n — m > 0, 2k — n > and 2j — m > 0. We conclude that 
small k and j in turn imply small n and m. The meaning of this will become clear after 
developing explicitly the first orders in the temperature expansion. 

To start with, let us focus on the simplest case, assuming that the temperature 
is low enough so that we need to keep only the (ksT) order, that is k — 0. This in 
turn implies n = m = 0. Then the suitable deformations are limited to those that are 
independent of g. We get 

6Q = jdsds'Ms)\ {j^{s)5{s-s')+f w {s,s')^%{s') (23) 

where we have used F° g = 1. We see that we have recovered the expression ([7]) for the 
zero temperature case. 

The next to leading order in temperature, (k B T) 2 , gives 

r 2 

5n = jdsds' £ i> n {s)M nm {s,s')i, m {s') (24) 

n,m=0 
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where M nm (s,s') is a 3 x 3 symmetric matrix given in terms of the J~ 1 (s)'s and the 
fnm(s, s')'s. M can be easily diagonalized to give 

5Q = jdsds' (Ms)M 00 Ms') + Ms)M n Ms')) (25) 

where i[>o(s) and ipi{ s ) are new arbitrary functions (which are linear combinations of the 
original ipo{ s ) an d ipi(s)), and the diagonal matrix elements read 

Moo = ~ (V(*) + AkT) 2 ^(a)) 5( S - s') + i/ o(s, s') + ^(kT) 2 f 02 (s, s') 

M u = ^p^J \ S )6(s-s') (26) 

where we have consistently expanded the matrix elements to order {kT) 2 , an we omitted 
the third eigenvalue M 33 since it scales as kT A . Replacing into the energy we get 

Sn^dsds'^^Jo^WikTfj^is)) S(s-s') + f o(s,s') + ^-{kT) 2 fo 2 {s, s')j Ms') + 

+ fi^l. Jdsds'Ms) J-\ s )5(s-s') Ms') (27) 

Note that the second line containing the ip\(s) modes, is proportional to the energy 
of a free fermion at zero temperature with a positive proportionality constant. This 
observation makes evident that no instability can be originated from those modes, and 
that we can safely turn them off. Focusing on the first line that contains the ^o( s ) 
modes, we see that the expression is analogous to ([7]) but with modified Jacobian and 
interaction function. Then we can directly apply the method of section [2] to analyze the 
instabilities. 

Corrections to the quadratic form coming from the higher orders in the temperature 
expansion can be obtained systematically in a similar way. 

4. Example: s-wave interaction in a square lattice at finite temperature 

In this section we apply the formalism just developed to a model of fermions in a square 
lattice with a s-wave interaction, whose zero temperature case was studied before in 
[50] . The comparison with the finite temperature case is interesting by itself, but the 
main aim of the present section is to show an example of application of the general 
procedure presented above. 

We start by considering free fermions with dispersion relation given by 

e(k) = -2t(cos(k x ) + cos(k y )) (28) 
At T = the FS is defined by 

g(k) =n- e(k) = (29) 

Now we change variables according to (T5]) in the following way 

g(k) =fj, + 2t(cos(k x ) + cos(k y )) (30) 
/ tan(k y /2) \ 



s(k) = arctan — ^ 

v ; \tan(k Q 



/2), 



(31) 
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• kT = 0.01 • kT = 0.5 kT = 1 

Figure 1. (Color online) a) Instability regions in the (//, Uq) plane for different values 
of kT. The dotted line represents the T — limit 30 . At T = the system becomes 
more unstable than at finite temperature, as it should be expected. The shaded valleys 
around /i — are outside the region of validity of our low T approximation, b) 
Instability regions in the (fx, kT) plane, the boundary of the colored region defines 
a critical value T e (/z) where the instability is reached. The shaded triangular region 
around (J, = above the main diagonal of the \x-kT plane, does not satisfy our low T 
criterion. 



The Jacobian of this transformation is given by 



J-\g, s) = 2ty/l-P(g-fj,) cos 2 (2s) (32) 

where 



x x 2 



fi{x) = 1 - (^-J (33) 
This Jacobian can be expanded in powers of g as in eq. (ITS]) , obtaining 

J~ l (s) = 2tyJl-P(fi) cos 2 (2s) (34) 

Jr\s) = - (35) 



4Jn^ 1 ( , 

(36) 



! _ cos 2 (2 g ) / _ /iW(2£)\ 
2 [ ) " ±Ja\s) V 4J - l2 ( S ) J 
In the case of the s-wave Pomeranchuk instability studied in [30], the interaction 
function for the Fermi liquid is given by the expression 

f(k,k') = U (37) 

Uq being a constant that measures the strength of the interaction. 

This form of the dispersion relation and interaction function can be obtained by a 
mean field approximation or a first order perturbative expansion of the Hubbard model 
in the square lattice, considering only hopping to nearest neighbors [24"1 125] 

^ = -*EE 4 A* + (2^) 2 ^o E (38) 

where (ij) indicates a sum over nearest neighbors. 



Pomeranchuk instabilities in lattice systems 



10 



We will expand to quadratic order in temperatures, so that the bilinear form we 
have to study is given in expression fl25l) . that specialized for the interaction f )37|) reads 

5Q = ~ J dsds' Ms) ((^(s) + 7r 2 (kT) 2 J^(s)) 5{s - s') + f/ ) Ms') (39) 

It can be shown that higher (quartic) orders in temperature correcting the 5(s — s') 
coefficient in the above expression can be safely discarded as long as kT <C //, while the 
corresponding corrections proportional to the interaction strength Uq are irrelevant as 
long as kT <C t. 

As explained above, we now choose an arbitrary basis of L 2 [FS], that in the present 
case will be the set of trigonometric functions {sin(ns), cos(ns)}. The convenience of 
this choice can be seen by the fact that each of the coefficients in the expansion (118p 
can be written as 

J-\s)=Y J 3 i n ) ^{Aks). (40) 

k 

Note that, if the cos(Os) = 1 mode has to be represented as a total derivative with 
respect to s, then 1 = d s X implies A = s which, being multivalued, does not belong to 
L 2 [FS]. This in turn implies that we are including a mode that does not satisfy the 
condition imposed by the Luttinger constraint. The deformation of the FS originated 
in that mode corresponds to the addition of particles to the system. 

Next we go through a Gram-Schimidt orthogonalization procedure to obtain an 
orthogonal basis {£i( s )}> whose pseudo-norms we evaluate. If for some i we 

get a negative pseudonorm, a deformation parameterized by £j(s) will have negative 
energy, and we can conclude that we have detected a Pomeranchuk instability in the 
i-th channel. 

We have studied the first 40 channels detecting which ones are unstable, and 
we plotted the resulting phase diagram in the (/i,kT,U ) space in Fig.l. There the 
instability region determined by the dominant unstable channels is displayed. 

When the temperature is increased at a constant value of the interaction Uq the 
system becomes more stable. This behavior is sketched in Fig. la where the (//, Uq) plane 
is shown for different values of kT, the colored areas corresponding to unstable regions. 
The phase diagram Uq vs. \i of [30] is recovered at T = 0. 

In Fig.l we show the Uq = —0.25 plane. The boundary of the instability region 
is strongly dependent on the chemical potential [i and we can infer the presence of an 
optimal value \x c where the Fermi liquid becomes unstable at the highest temperature. 

In Figj2] we show a tree dimensional phase diagram in the space (fi,kT,Uo), in 
which the above described features can be seen in detail. 

The interesting features of the phase diagram obtained by applying the method to 
the present model, makes evident its potential power to study instabilities of the ci-wave 
type or more complex interactions, as those discussed in [22J in the context of URu2Si2, 
etc. 

We have analyzed only the first 40 modes for illustration purposes. The scaling of 
the unstable region as a function of the mode index i can be studied as in [30] and [31] 
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Figure 2. (Color online) Instability region in the (fi,Uo,kT) space. The shaded 
central region around the /i = plane is outside the range of validity of our low T 
approximation. 



to show that the region of instability is bounded. 

5. Extension of the method to systems with spin/color 

In the present section we will enlarge the approach to include the case of multiple 
fermionic species (like spinful fermions, multiband materials, etc.). Since in such 
materials we have one FS for each species, we will denote them with a new (internal) 
index a, running from 1 to N (e.g. in the spin 1/2 case, N = 2). Then the equations 
defining the FS's are given by 

/2 a -e a (k)=0 a = l,...,N (41) 

Now the constraint imposed by the Luttinger theorem on Sn a (k) reads 

N 

/ d\8n a (k) = 0. (42) 

The changes of variables have to be done on each FS, to go from (k x ,k y ) to normal 
Fermi variables (g a ,s a ): 

g a = g a (k) = ft a -e Q (k), 

s a = s a {k) , (43) 
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where g a is normal to the Fermi Surface corresponding to species a and the tangent 
variable s a runs between — 7r and ir. The associated Jacobians are given by 

d{g a , s° 



(44) 



dk 

and then the constraint can be rewritten as 

N 

/ dg a ds a J a (g a , s a )5n a (g a , s a ) = . (45) 

Now writing as before 

5n a (g, s) = H[g a + 5g a (g, a)] - H[g a ] , (46) 
and replacing it into the constraint (115]) we get 

N r r s 9 a (9 a ' sa ) 

£ ds a dg a J a (g a ,s a )=0. (47) 

a=l J J ° 

Performing the integrals over the g a, s to lowest order in 5g a 

N 

/ ds a J a (s a )5g a (s a ) = . (48) 

where J a (s a ) = J a (g a , s a )\ ga=Q and 5g a (s a ) = 5g a (g a , s a )\ g c =0 . Being the s a 's dumb 
variables, we can use the same name on each FS and rewrite the previous expression as 
a single integral over one variable s 

J ds (y2j a (s)Sg a {s)\ =0. (49) 

which can be solved as before by writing 

N 

J2 Us)5g a (s) = d s \(s) (50) 

a=l 

with A(s) an arbitrary function of s. Up to this point, we have merely added spin/color 
indices to the procedure outlined in Section |2j but now we note that the above expression 
can be further simplified by writing 

Wis) = tL (d s X(s)J a (s) + Ji{sj) (51) 

where |J| 2 = J2 a Ja{s) J a (s) , and J„(s) is an arbitrary vector in the internal space 
of a indices that satisfy J2 a J a ( s )Ja(s) = 0. This J„(s) represents the new degrees 
of freedom that were absent in the spinless case, that correspond to excitations that 
change the spin/color of the quasiparticles. We will need N — 1 new arbitrary functions 
to parameterize them. 

In the particular case of spin 1/2, in which we will concentrate in what follows, we 
have a = 1, 2 corresponding to the spin up and spin down FS's respectively Any vector 
perpendicular to J a can be written as 

Ji = v{s)E^ J P> (52) 
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where e a/3 is the Levi-Civita antisymmetric symbol, and r](s) is the new arbitrary 
function of s needed to parameterize the part of 5g a perpendicular to J a . 

Let us now compute the variation of the ground state energy corresponding to 
variations on the occupation numbers 5n a (k) 

N N 1 , , 

5Q=J2 d\{e a {k)-fi a )5n a {k)+ ]T - d\ d\'f a ^(k,k') Sn a (k)8n^(k') . (53) 
We now change variables as in (|43|) and express 5n a (k) using (|46l) to get 

N ,0 , 

6Q = -Y, d 9 a ds a Ug a ,s a )g a 

~, J-5a a J 



a=l ■ 



N 



+ 1 E 1° d 9 a f dg^d S a [d S ^r^g a }S a } g^ S ^J a (g a ,s a )Mg a }S a ). 

(54) 

Expanding to lowest order in Sg a and renaming the dumb variables s a we get 

r r N 1 
5Q = JdsJds'Y, Sg a (s)5g^s')-(j a (s)6 a ^(s - s') + f a ^s,s') J a (s)J p (s')) (55) 

a, /3=1 ~ 

where f a/3 (s,s') = f a/3 (g a ,s,g'P,s')\ ga=g ,f3 =0 . Other than the additional color/spin 
indices, the crucial difference of the above result from the formula obtained in section 
|2] will become evident after the explicit solution for the Sg a, s which are allowed by the 
Luttinger theorem (IST!) are inserted. We then obtain 

5Vl= fdsfds' Yl ii a (s)K ab (s,s')i}\s'). (56) 

J a,b=l,2 

Here the new indices a, b = 1, 2 are not spin/color indices, but refer to the two arbitrary 
functions that parameterize the FS deformations (as in floTl) and fl52l) ). namely 

4>\s) = d s X(s) , ^ 2 ( s )= v (s) (57) 

and the kernel K (s, s') is consequently a two-by-two matrix, which entries are given 
explicitly by 

K "(*> s ') = oi n.w THi^ IE J l^ s ~ s ') + E Ja(s)f aP (s, s')4(s') 
R22 ^ s ') = oi r ^M2i 7^12 ( E ^^JaWi* - s') + 



2\J(s] 


>IW)| 2 \ 




>IW)I 2 ( 


2\J(s] 




+ 













2\J(s] 


>IW)| 2 ( 






2\J(s] 


>IW)| 2 \ 



K 12 (s, s') = ol7 ^, 2 , 7 ^ AI2 ( E ^JpS(s - s') + £ e*JZ(8) J,(s')J s (s')n(s, s') 
K ^ <0 = 777777777777^ ( E ~ + 
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+ £ J a {s)J p {s) J 5 (s')J s (s')f aS (s, s')) (58) 

This is the kernel to be diagonalized by the Gram-Schmidt procedure. The main 
difference then is that we get an extra matrix structure in the two dimensional space of 
deformations, parallel and perpendicular to the Jacobian J a , and labelled by the indices 
a, b. This corresponds to separating those deformations that preserve total spin/color 
from those which do not. 

The extension just presented of the generalized Pomeranchuk method can be used 
to study the instabilities of the Fermi liquid under perturbations that involve spin flips, 
or induced by the presence of an external magnetic field. 

6. Summary and Conclusions 

In the present work we have generalized the extended Pomeranchuk method developed 
in [30] and [31] to detect instabilities in more realistic lattice systems. We have extended 
the method along two main directions: 

First, in order to get closer to experiments, we have modified the method to include 
finite temperature effects. This allows to study the extended phase diagram and in 
particular, to compute the critical temperature as a function of the chemical potential 
and the interaction strength. 

As a second direction we modified the method in order to make it applicable to study 
multiband systems including spin or any other internal degree of freedom. We show 
that in those cases an extra matrix structure appears, parameterizing the deformations 
according to whether they preserve or not the associated quantum number. 

As an example of application, the effects of a finite temperature on the phase 
diagram of the fermion system in the square lattice with an s-wave instability were 
analyzed. This calculation shows the simplicity of the method and its potentiality to 
identify the unstable regions of the phase diagram. 

After these generalizations, the method is now ready to be applied to models 
relevant for different novel materials which are thought to realize different types of 
Pomeranchuk transitions, such as the bilayer ruthenate S^R^Oy, the hidden order 
transition material URu2Si2, etc. We expect to study these problems in the near future. 
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